arXiv: 1504.03172v2 [physics.bio-ph] 16 Apr 2015 


Direct measurement of DNA-mediated adhesion between lipid bilayers 


S. F. Shimobayashi 

Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan 

B. M. Mognetti 

Interdisciplinary Center for Nonlinear Phenomena and Complex Systems & Service de 
Physique des Systemes Complexes et Mecanique Statistique, Universite libre de Bruxelles, 
Campus Plaine, CP 231, Blvd du Triomphe, B-1050 Brussels, Belgium. 

L. Parolini, P. Cicuta,* and L. Di Michele^ 

Biological and Soft Systems Sector, Cavendish Laboratory, 

University of Cambridge, CB3 OHE, Cambridge, United Kingdom 

D. Orsi 

Department of Physics and Earth Sciences, University of Parma, 1-43124, Parma, Italy 

(Dated: April 20, 2015) 

Multivalent interactions between deformable mesoscopic units are ubiquitous in biology, where 
membrane macromolecules mediate the interactions between neighbouring living cells and between 
cells and solid substrates. Lately, analogous artihcial materials have been synthesised by functional¬ 
ising the outer surface of compliant Brownian units, for example emulsion droplets and lipid vesicles, 
with selective linkers, in particular short DNA sequences. This development extended the range of 
applicability of DNA as a selective glue, originally applied to solid nano and colloidal particles. On 
very deformable lipid vesicles, the coupling between statistical effects of multivalent interactions and 
mechanical deformation of the membranes gives rise to complex emergent behaviours, as we recently 
contributed to demonstrate [Parolini et al. Nature Communications, 2015, 6, 5948]. Several aspects 
of the complex phenomenology observed in these systems still lack a quantitative experimental 
characterisation and fundamental understanding. Here we focus on the DNA-mediated multivalent 
interactions of a single liposome adhering to a flat supported bilayer. This simplified geometry en¬ 
ables the estimate of the membrane tension induced by the DNA-mediated adhesive forces acting on 
the liposome. Our experimental investigation is completed by morphological measurements and the 
characterisation of the DNA-melting transition, probed by in-situ Forster Resonant Energy Transfer 
spectroscopy. Experimental results are compared with the predictions of an analytical theory that 
couples the deformation of the vesicle to a full description of the statistical mechanics of mobile 
linkers. With at most one fitting parameter, our theory is capable of semi-quantitatively matching 
experimental data, confirming the quality of the underlying assumptions. 


Due to both fundamental and practical interest, the 
complex phenomenology of multivalent selective inter¬ 
actions has recently been a hot research topic in Soft 
Matter and Physical Chemistry. One of the main driving 
forces behind this effort has been the development of 
self-assembly strategies based on DNA-mediated multi¬ 
valent interactions. [IJIntroduced by the seminal works of 
Mirkin[2] and Alivisatos,[3] who nearly two decades ago 
demonstrated the selective DNA-mediated self-assembly 
of gold nanoparticles, this approach has been optimised 
to the point of mastering the structure of multicompo¬ 
nent crystal lattices[4-8] and amorphous materials.[9-11] 
Applications of these self-assembly strategies span from 
photonics and plasmonics,[12] to biosensing[13], and gene 
therapy. [14-16] An abundance of experimental studies 
called for the development of models to understand 
the complex statistical mechanics of DNA-mediated 
multivalent interactions between solid particles. [17-22] 
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More recently, the same self-assembly strategy has been 
applied to compliant Brownian units, including emulsion 
droplets, [23, 24] and lipid vesicles. [25-29] Especially 
for the case of liposomes, which are significantly more 
deformable than emulsion droplets, DNA-mediated 
adhesion plays at a similar magnitude as other forces, 
leading to rich coupling. [28, 30] Moreover, on these 
liquid surfaces, DNA linkers can freely diffuse, and the 
confinement induced by binding results in significant 
entropic contributions to the hybridisation free energy. 
We recently demonstrated that the coupling between 
mechanical deformability of Giant-Unilamellar-Veisicles 
(GUVs) and DNA-mediated interactions, results in 
unexpected emergent response to temperature changes, 
leading to negative thermal expansion and tuneable 
porosity of DNA-GUV assemblies. [30] By means of an 
analytical theory, that jointly describes the geometrical 
features of the GUVs and the statistical mechanics 
of DNA tethers, [30, 31] we rationalised experimental 
evidence and highlighted the importance of translational 
entropy in coupling the DNA-binding free energy to the 
morphology of the substrates. Nonetheless, a further 
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experimental investigation is necessary to assess all the 
aspects of the complex phenomenology of DNA-GUV 
conjugates, and test the accuracy of the theoretical 
framework applied to their description. In particular, 
a direct measurement of the DNA-mediated adhesive 
forces and the quantitative characterisation of the 
melting transition would represent new strong tests of 
the current understanding. 

In this article, we present experiments aimed at mea¬ 
suring the temperature-dependence of the membrane 
tension induced on giant liposomes by DNA-mediated 
adhesion. Membrane tension is accessed by flickering 
measurements, in which we reconstruct the spectrum of 
thermal fluctuations of the liposomes at the equatorial 
cross-sections. [32-34] This is a classic technique widely 
used to study vesicle tension and bending rigidity. To 
make these measurements possible, we adopt an exper¬ 
imental geometry in which liposomes do not interact 
with each other, instead they adhere to Supported Lipid 
Bilayers (SBLs) fabricated on rigid glass substrates, 
as demonstrated in Fig. 1. This allows imaging of the 
equator of the vesicle through confocal microscopy. 
Various morphological observables are also precisely 
assessed via 3-dimensional confocal imaging. Finally, 
the use of suitably labelled DNA tethers enables the 
direct assessment of the relative number of formed DNA 
bonds by means of in-situ Forster Resonant Energy 
Transfer (FRET) spectroscopy. Experimental results 
are compared to predictions from our analytical theory, 
extended from ref. [30] to apply to this new geometry. 
Using at most one fitting parameter, we find semi- 
quantitative agreement for all the measured quantities, 
confirming the accuracy of the underlying assumptions 
of our model. 

The remainder of this article is structured as follows. 
First we describe the experimental setup, including 
sample preparations protocols, materials, imaging and 
image analysis techniques. Then we outline our theoret¬ 
ical model focusing on the novel aspects introduced here 
for the GUV-plane geometry. Finally we discuss the 
experimental results and compare them with theoretical 
predictions. 


I. EXPERIMENTAL 

Figure 1 shows a schematic of a DNA-functionalised 
GUV adhering to a DNA-functionalised SBL. GUVs and 
SBLs are prepared and separately functionalised with 
cholesterol labeled DNA constructs. The hydrophobic 
cholesterol inserts into the lipid bilayer allowing DNA 
tethers to freely diffuse. Gonnected to the cholesterol 
anchor there is a section of length L = 9.8 nm (29 base 
pairs[35]) of double-strands DNA (dsDNA), terminating 
in a single-stranded DNA (ssDNA) sticky end, which me¬ 
diates the attractive interactions. To further facilitate 
the pivoting motion of the DNA tethers, 4 unpaired ade- 



FIG. 1. Not-to-scale schematics of a DNA-functionalised 
GUV adhering to a supported bilayer and details of the func¬ 
tionalisation. 


nine bases are left between the cholesterol anchor and 
the dsDNA spacer. One unpaired adenine is left be¬ 
tween the dsDNA spacer and the sticky end. Two mu¬ 
tually complementary sticky ends are used: a and a', of 
7 bases each (see Fig. 1). GUVs and SBLs are function¬ 
alised with equal molar fractions of both a and a'. Teth¬ 
ers can therefore form intra-membrane loops and inter¬ 
membrane bridges. The latter are responsible for the ob¬ 
served adhesion and are confined within the contact area 
between the GUVs and the SBLs. We use FRET spec¬ 
troscopy to estimate the fraction of formed DNA bonds. 
To enable these measurements the termini of the sticky 
ends a and a' are functionalised with a Gy3 (indocarbo- 
cyanine) and a Gy5 (indodicarbocyanine) molecules re¬ 
spectively, spaced by an unpaired adenine base. Note 
that the sticky ends used here are shorter than those 
adopted in our previous study[30] (7 base-pairs instead 
of 9), this choice was made to lower the melting temper¬ 
ature of the DNA to within an experimentally accessi¬ 
ble temperature, and reduce the overall strength of the 
DNA-mediated adhesion, making it easier to measure by 
flickering spectroscopy. 


A. Materials and sample preparation 

GUVs electroformation 

l,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPG, 
Avanti Polar Lipids) GUVs are prepared by standard 
electroformation in 300 mM sucrose (Sigma Aldrich) 
solution in double-distilled water. [36, 37] Full details can 
be found in ref. [30] 
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1. Supported bilayers 

DOPC supported bilayers (SBLs) are formed by rup¬ 
ture of Small Unilamellar Vesicles (SUVs ~ 100 nm) on 
the hydrophilised glass bottom of sample chambers. [38, 
39] SUVs are prepared by standard extrusion. Briefly, 
200 /il of 25 mg/ml DOPC solution in chloroform are 
left to dry in a glass vial, then hydrated by adding 500 
fA of 300 mM sucrose solution and mixed by vortexing 
for at least 5 minutes. The solutions are then trans¬ 
ferred in plastic vials and treated with 5 cycles of rapid 
freezing/unfreezing by alternatively immersing the vial in 
baths of liquid nitrogen and warm water. [27] Extrusion 
is carried out using a hand-driven mini-extruder (Avanti 
Polar Lipids) with a polycarbonate track-etched mem¬ 
brane (100 nm pores, Whatman). To facilitate rupture 
of the SUVs and bilayer formation, the extruded solution 
is diluted in a 1:9 ratio in iso-osmolar solution containing 
TE buffer (10 mM tris(hydroxymethyl) aminomethane, 
1 mM ethylenediaminetetra acetic acid. Sigma Aldrich), 
5 mM MgCl 2 and 272 mM glucose (Sigma Aldrich). 
Sample chambers are obtained by applying adhesive 
silicone-rubber multi-well plates (6.5 mm x 6.5 mm x 3.2 
mm, Flexwell, Grace Biolabs) on glass coverslips (26x 60 
mm, No.l, Menzel-Glsaer), cleaned following a previously 
reported protocol. [9] To form SBLs, the glass bottom of 
the cells is hydrophilised by plasma cleaning on a plasmo- 
chemical reactor (Femto, Diener electronic, Germany), 
operated at frequency of 40 kHz, pressure of 30 Pa, and 
power input of 100 W for 5 minutes. Within 5 minutes 
from plasma cleaning, each cell is filled with 100 pil of 
diluted SUV solution and incubated for at least 30 min¬ 
utes at room temperature to allow bilayer formation. To 
remove excess lipid and magnesium, the cells are repeat¬ 
edly rinsed with the experimental solution (TE buffer, 
100 mM NaGl, 87 mM glucose). Gare is taken to always 
keep the bilayers covered in buffer to avoid exposure to 
air. 

For experiments not involving FRET spectroscopy 
GUVs and SBLs are stained with 0.8-1.2% molar 
fraction of Texas Red l,2-Dihexadecanoyl-sn-Glycero-3- 
Phosphoethanolamine, Triethylammonium Salt (Texas 
Red DHPE, Life Technologies). 

DNA preparation 

The DNA tethers are pre-assembled from two ssDNA 
strands, one of them (i) functionalized with a choles¬ 
terol molecule and the second (ii) carrying the sticky end: 

i 5'-CGT GCG CTG GCG TCT GAA AGT 

CGA TTG CG AAAA [Gholesterol TEG] -3' 

it 5-GG GAA TCG ACT TTG AGA CGC 

GAG CGC ACG A [Sticky End] A Gy3/Gy5 -3'. 

The bold letters indicate the segments forming the ds- 
DNA spacer, the italic letters the inert flexible spac¬ 


ers. DNA is purchased lyophilized from Integrated DNA 
Technologies, reconstituted in TE buffer, aliquoted, and 
stored at -20 °G. For assembling the constructs, we di¬ 
lute equal amounts of the two single strands, i and ii, 
to 1.6 /iM in TE buffer containing 100 mM NaGl. Hy¬ 
bridization is carried out by ramping down the temper¬ 
ature from 90°G to 20°G at a rate of -0.2°G min“^ on a 
PGR machine (Eppendorf Mastercycler). [30] 


Membrane functionalisation. 


Fuctionalisation of the supported bilayers is carried 
out by injecting 90 ptL of iso-osmolar experimental 
solution (TE buffer, 87 mM glucose, 100 mM NaGl) 
containing overall X moles of DNA constructs into each 
of the silicone-rubber cells, with equal molarity of a 
and a' strands. Similarly, GUVs are functionalised by 
diluting 10 ^L of electroformed vesicle solution in 90 
/iL of iso-osmolar experimental solution containing X 
moles of DNA constructs. GUVs and SLBs are incu¬ 
bated for at least 1 hour at room temperature to allow 
grafting. After incubation, 10 /il of the the liposome 
solution are injected into the sample chambers, which 
are immediately closed with a second clean coverslip and 
sealed using rapid epoxy glue (Araldite). Gare is taken 
to prevent the formation of air bubbles. Sedimentation 
of the GUVs results in the formation of an adhesion 
patch between GUVs and supported bilayer. In a typical 
sample a fraction of the GUVs is found to form clusters. 
We limit our analysis to isolated GUVs. 

We tested samples at different DNA concentrations, 
obtained by setting X to Aiow = 0.09, Amed = 0.9 and 
Ahigh = 2.7 pmoles. We have previously quantified the 
surface coverage of GUVs functionalised with Amed = 0.9 
pmoles in p^NA = 390 ± 90 /rm“^.[30] Here we propor¬ 
tionally assume = 39±9 /tm“^, pJ],na = 1200 ±300 
pm“^. Having estimated that the overall surface of 
the SBL is approximately equal to the surface of the 
GUVs used for each sample, we assume equal DNA 
coverages for the SBL. This assumption is confirmed by 
fluorescence emission measurements: at sufficiently high 
temperature, when no DNA bridges are formed between 
GUVs and SBL, DNA is uniformly distributed on both 
interfaces. In this regime the fluorescence emission from 
DNA located within the contact area between a GUV 
and the SBL approximately equals twice the intensity 
measured on the free SBL (2.1 ± 0.1), confirming that 
GUVs and SBL have, within experimental errors, the 
same DNA coverage. 
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B. Imaging and image analysis 


1. Flickering analysis 


Imaging and temperature cycling. 


The samples are imaged on a Leica TCS SP5 laser¬ 
scanning confocal microscope. Texas Red DHPE is 
excited with a He-Ne laser (594 nm). For FRET spec¬ 
troscopy measurements, Cy3 is excited with an Ar-ion 
laser line (514 nm). The temperature of the sample is 
regulated with a home-made Peltier device controlled 
by a PID (proportional-integral-derivative) controller, 
featuring a copper plate to which the sample chamber is 
kept in thermal contact. Two thermocouples are used 
as temperature sensors. The first sensor, kept in contact 
with the copper plate, serves as a feedback probe for the 
PID controller. The second thermocouple is inserted in 
a dummy experimental chamber, filled with water, and 
used to precisely probe the temperature of the sample. 
For all the temperature-dependent experiments imaging 
is carried out using a Leica HCX PL APO CS 40x 0.85 
NA dry objective, to prevent heat dissipation. 

The temperature-dependent morphology of adhering 
GUVs is captured via confocal z-stacks and reconstructed 
using a custom script written in Matlab. Briefly: Each 
z-stack contains a single adhering GUV. The plane of 
the SBL/adhesion patch is identified as the one with 
maximum average intensity, then the adhesion patch is 
reconstructed by thresholding, following the application 
of a bandpass filter to flatten the background and remove 
pixel-level noise. From the area Ap of the adhesion 
patch we extract the patch-radius as Rp = Apj-n. The 
portion of the z-stack above the SBL is then scanned, 
and in each plane the contour of the vesicle, suitably 
highlighted by filtering and thresholding, is fitted with a 
circle. The slice featuring the largest circle is identihed 
as the equatorial plane, determining the vesicle radius 
R. The contact angle is derived 6 — sin~^{Rp/R) 
(see Fig. 1). 

For flickering experiments, movies are recorded across 
the equatorial plane. Details of the flickering analysis 
are reported in the next section. To improve the quality 
of the signal, imaging for morphological characterisation 
and flickering analysis is carried out on samples stained 
with Texas Red DHPE. 

FRET spectroscopy measurements are carried out 
by performing a spectral emission scan (A-scan) of 
the contact area between an adhering GUV and the 
surrounding free SBL. While exciting the donor (Cy3), 
the emission of donor and acceptor is reconstructed 
by scanning the acquisition window from 530 to 785 
nm, with a 6.75 nm binning. Similarly to the case of 
z-stacks, the adhesion patch is identified in each image 
by filtering and thresholding. For FRET imaging, non- 
fluorescent lipids are used in order to prevent undesired 
energy transfer between the bilayer and the DNA tethers. 



FIG. 2. Intermediate steps of the contour detection algorithm, 
a. The position of the centra and the mean radius R of the 
vesicle define an annular region in the frame; b, the annular 
region C{r, ip) extrapolated from the frame is used to compute 
an averaged template {C{r)}; a 2D correlation Q{r,<j)) of the 
two is computed and fitted to find the membrane position 
with sub-pixel precision, as described in the text (steps 1-6). 
The resulting contour is reported as a black line; c, colormap 
representation of the membrane fluctuations relative to the 
average radius. The colorbar reports the relative fluctuation 
as percentage in the range —5%, -|-5% 

Typically, in case of phase-contrast imaging, the con¬ 
tour of fluctuating GUVs is reconstructed by finding 
the inflection point in the radial intensity profile, as in 
reference. [34] In case of fluorescence imaging, the maxi¬ 
mum of the intensity proHle is commonly used to mark 
the membrane position. Here, we designed a Matlab al¬ 
gorithm to reconstruct the contour of GUVs from their 
equatorial cross section with sub-pixel precision, indepen¬ 
dently on the imaging mode employed. The algorithm 
proceeds following these steps: 

1. The position of the centre, and a rough-guess value 
for the radius R of the GUV are manually selected 
by the user on the first frame of the video (see 
Fig. 2a). 
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2. The image of the radial profile C{r, cp) is selected 
within an annular region of user-defined width that 
contains the membrane. The annular region is then 
mapped onto a rectangular stripe using a cubic in¬ 
terpolation (see Fig. 2b, top). 

3. A template (C'(r)) of the radial intensity profile is 
calculated from the annular region, either by man¬ 
ual selection, or by automatically averaging along 
(p (see Fig. 2b, top-right). 

4. The algorithm computes the 2D correlation func¬ 
tion Q{r,(p) of C{r,(p) with the template {C{r)). 
For each value of the polar angle tp, Q(r, tp) is max¬ 
imum at the radial coordinate r where the tem¬ 
plate best correlates with the membrane profile (see 
Fig. 2b, bottom). 

5. A rough measurement of the contour position r{ip) 
is obtained by fitting the peak of the correlation 
image with a parabola around its maximum, for 
every value of (p. The algorithm cycles three times 
from point 2 to 5 , each time refining the estimate 
of the vesicle’s centre and mean radius using the 
obtained contour. 


6. The rough contour estimate is used as the initial 
point of a global minimisation of the energy-like 
functional 


e(A V?) = XI 


[1 -Q{r{ip),(p)]+a 


d‘^r{(p) 


( 1 ) 


here, the membrane position corresponds to the 
minimum of 1 — Q, corrected by a term that gives 
low weight to high curvature features, commonly 
due to random noise in the image. This retrieves 
the measurement of the contour position r{(p). 
The minimisation makes use of the Nelder-Mead 
algorithm [40]. 


The procedure is repeated for every frame in a video of 
the fluctuating membrane, each time using the center and 
radius values of the previous frame as the starting point 
of the algorithm. An example of the temporal evolution 
of the membrane radial profile r{ip) is shown in Fig. 2c. 
The spectrum of the thermal fluctuations of the contour 
is calculated using Matlab Fast Fourier Transform algo¬ 
rithm, using the formula 




where M is the number points used to map the contour, 
and then averaged over the ensemble of frames. 


2. Control experiments 

Control experiments are performed to measure 
temperature-dependent membrane tension of non¬ 
adhering GUVs. Plain, non-functionalised, GUVs are 


imaged while laying on a glass substrate passivated with 
bovine serum albumin (BSA, Sigma Aldrich). Part of 
the images is carried out in bright-held microscopy, us¬ 
ing a Nikon Eclipse Ti-E inverted microscope, a Nikon 
PLAN APO 40x 0.95 N.A. dry objective and a IIDC 
Point Grey Research Grasshopper-3 camera. For con¬ 
trol experiments GUVs are diluted in a 1:9 ratio in iso- 
osmolar glucose solution to enable sedimentation. 


II. THEORETICAL MODEL 

We present a quantitative model describing the DNA- 
mediated adhesion of GUVs on flat supported bilay¬ 
ers. The model is adapted from reference, [30] where we 
treated the case of two identical adhering GUVs. Let us 
consider the interaction between an inhnite SBL and a 
GUV adhering to it 

U = Umembrane T UdNA T Uq, (2) 

where Umembrane accounts for the mechanical deforma¬ 
tion of the GUV, Udna encodes for DNA-mediated ad¬ 
hesion, and Uq is the reference energy, calculated for iso¬ 
lated GUV and SBL. 


A. Membrane deformation 


In Eq. 2, Umembrane Summarises three contributions: 
stretching energy, bending energy, and the entropic cost 
of suppressing thermal fluctuations of the contact area 
between the GUV and the SBL. [41] In the limit of strong 
adhesion the stretching energy dominates over the other 
two contributions, which we can neglect.[41] As discussed 
later, for the case of DNA-mediated interactions, this 
condition is generally verified at low enough temperature. 
We can rewrite 


f^membrane (^) ^stretching (^) 


A{e) - A 


2A 


( 3 ) 


where ATa is the stretching modulus of the membrane, 
A(0) is the overall (stretched) area of the GUV, and A is 
the reference, un-stretched, area. In the limit of strong 
adhesion the GUV will take the shape of a truncated 
sphere with contact angle 9 (Eig. 1), which we take as 
the independent variable of our model. [41] Within the as¬ 
sumption of constant inner volume V = AttHq/H, where 
Rq is a reference radius, the total area A of the GUV 
and the adhesion patch area Ap can be expressed as a 
function of the contact angle 9 (see Appendix, section 
Al). The un-stretched vesicle area A exhibits a strong 
temperature-dependence 

i = Ao [1 + a(T - To)], (4) 


where a is the area thermal expansion coefhcient [42] and 
To is the neutral temperature of the GUVs, at which its 
reduced volume equals unity and A = 47ri?Q (see Ap¬ 
pendix, section A 1). 
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B. DNA-mediated adhesion 


where 


We now focus on the DNA-mediated contribution to 
the interaction energy in Eq. 2: C/dna- 
Given that the persistence length of dsDNA is ~ 
50 nm ^ L = 9.8 nm,[35] we can model the dsDNA 
spacers as rigid rods that, thanks to the fluidity of 
the membrane, can freely diffuse on the surface of the 
bilayers. [30, 31] Free pivoting motion is guaranteed by 
the flexibility of the joint between the cholesterol anchor 
and the dsDNA spacer (Fig. 1). As demonstrated in 
ref. [30], we can regard the sticky ends as point-like re¬ 
active sites, neglecting their physical dimensions. More¬ 
over we can safely assume that the distance between the 
adhering membranes within the contact area is equal to 
L.[30] Finally, we neglect excluded volume interactions 
between unbound DNA tethers. [30, 43] The last two as¬ 
sumptions guarantee a uniform distribution of unbound 
DNA tethers and loops over the GUV and SBL surfaces. 
The free energy change associated to the formation of a 
single bridge (b) or a loop {£) within the GUV is 

AGb/, = AG° - rAA^°/f, (5) 

where AG° = Ai4° — T AA^ is the hybridisation free en¬ 
ergy of untethered sticky ends, which can be calculated 
from the nearest-neighbour thermodynamic model, [44] 
eventually corrected to account for neighbouring non- 
hybridised bases. [45] In Eq. 5, the term —TAS']]™^ ac¬ 
counts for the confinement entropic loss taking place 
when tethered sticky ends hybridise,[18, 20, 46] and can 
be estimated as (see Appendix, section A 2) 


AG(6i) = AG° 


feTlog 


1 

PoLA{9) 


( 6 ) 


where po = IM is a reference concentration, and we high¬ 
lighted the coupling with the geometry of the GUV via 
the d-dependence. Note that, contrary to the case of two 
adhering GUVs,[30] here AA™''^ = AS'™"^ The local 
roughness of the membranes could also influence AG [47] 
- this effects will be studied elsewhere. 

We indicate with N the total number of a and a' tethers 
on the GUVs, and model the SBL as an infinite reservoir 
of tethers. A combinatorial calculation detailed in the 
Appendix (sections A 3 and A 4) allows the derivation of 
the overall hybridisation energy for the system of linkers 


Uhyh = Nk^T 


X( -k 21og(l - xe - Xh) - 2 


Ni' 

Iv 


( 7 ) 


In Eq. 7, Xh and xi are the fraction of tethers involved 
in bridges/loops, given by 


_ 

{I -XI- a;b) 


( 8 ) 


XI 

{I -XI- Xhf 


( 9 ) 


qe = exp [-/3AG*], (10) 

( 7 b = ^ exp [-/3AG*]. (11) 

In Eqs. 10 and 11 the hybridisation free energy is re¬ 
defined as AG* = AG — kBTlogN. The combina¬ 
torial contribution —kBTlogN, typically estimated in 
~ —lOfceT, has a stabilising effect. 

The quantity Nf appearing in Eq. 7 indicates the average 
number of unbound DNA tethers anchored to the SBL 
available within the adhesion patch. The concentration 
Cf of unbound tetheres of each type (a or a'), such that 
N{{0) = CfAp(6<), is given by 


Co - Cf = Cf 


exp[-/3AGo] 

PqL 


( 12 ) 


where Cq = Pdna/2 is the total concentration of a and a' 
tethers. The full derivation of Eq. 12 is provided in the 
Appendix (section A 4). Note that the concentration of 
loops on the SBL is c^ = cq — Cf. 

Equation 7 generalises the expression found in ref. [30] to 
the case in which vesicles are in contact with an infinite 
reservoir of thethers. The overall DNA contribution to 
the interaction energy in Eq. 2 is 


C7dNA = Uhyh — 2NkBT\og 



(13) 


where the term on the right-hand side accounts for 
the change in overall confinement entropy following the 
area-change of the GUV. [30] 


C. Overall interaction energy. 

By combining Eqs. 3, 7, and 13 into Eqs. 2, we obtain 
an analytical expression for the interaction energy of the 
GUV-kSBL system as a function of the only independent 
variable 9 


U{9) -Uo = K^T 


A{9) - A 


2A 


NPbT 


Xi{9) + 


2 log [1 - Xi{9) - Xb(6')] - 2^^^ - 2 log 


Ao Jy 


(14) 


Note that U depends on 9 though the adhesion area 
Ap (see Eq. A21) and the total area A of the GUV. On 
the other hand, the reference energy Uq is ^-independent 
(see derivation in the Appendix, section A 5). The in¬ 
teraction energy can be minimised to calculate all the 
morphological observables (e.g. the contact angle 9, area 
of the adhesion patch Ap, area of the spherical section 
of the GUV A), as well as the fraction of formed bridges 
and loops (xb/^). 


= qe 
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FIG. 3. From left to right: 3D reconstruction from confo- 
cal z-stack, and confocal cross section of an adhering GUV 
(a) and a non-adhering GUV on a glass substrate (b). For 
3D rendering, confocal images have been acquired with a Le- 
ica HCX PL APO GS 63 x 1.4 NA oil immersion objective 
for better resolution, and deconvolved using the experimental 
point-spread function. Scale bars: 10 /rm. 



D. Model parameters and error propagation 


The model features seven input parameters: the ther¬ 
mal expansion coefficient a, the stretching modulus ATa, 
the length of the dsDNA tether L, the hybridisation en¬ 
thalpy and entropy AS® of the sticky ends, the 

DNA coating density pdna (used to calculate the overall 
number of strands per GUV - 27V), the neutral temper¬ 
ature To and radius Rq. 

The stretching modulus is estimated from literature data 
as ATa = 240±90 mN m“^, with the error bar covering the 
entire range of reported values. [48-51] The thermal ex¬ 
pansion coefficient has been experimentally estimated as 
a = 1.3±0.7 K“^.[30] The hybridisation enthalpy and en¬ 
tropy are estimated according to the nearest neighbours 
thermodynamic rules[44] as Ai7° = —54 ± 5 kcal mol“^ 
and AS^ = —154 ± 13 cal mol“^ K~^. It is not clear 
whether the stabilising effect of the non-hybridised dan¬ 
gling bases neighbouring the sticky ends, in the present 
case the adenine bases present at both sides of the se¬ 
quences (see Fig. 1), should be taken into account. It is 
indeed possible that their attractive contribution is com¬ 
pensated or overwhelmed by the repulsive effect of the 
long inert DNA connected to the sticky ends. [45] The er- 
rorbars in Ai7° and A5'° are included to cover both these 
scenarios. The length L = 9.8 nm of the dsDNA spac¬ 
ers can be precisely estimated from the contour length of 
dsDNA (0.388 nm per base-pair[35]). The DNA coating 
density is estimated for different samples as explained in 


the experimental section. The neutral temperature Tg 
changes widely from vesicle to vesicle due to the polydis- 
persity of electroformed samples. We generally use Tq as 
a fitting parameter. For a known Tq, the neutral radius 
Rq is experimentally determined as Rq = R\t= To ■ If Tq 
falls outside the experimentally accessible temperature 
range, we extrapolate 


Ro 


MTi) 


TT [1 + a (Ti - To)] ’ 


(15) 


where Ti is the minimum experimentally accessible tem¬ 
perature. 

Errors in the input parameters are numerically propa¬ 
gated to the theoretical predictions. [30] Briefly, we sam¬ 
ple the results of the model using random values of the 
input parameters X ± AX extracted from a Gaussian 
distribution with mean equal to X and standard devia¬ 
tion equal to AX. From a sample of size > 10000, we 
estimate the theoretical prediction of each observable as 
the median of the sampled distribution. Errobars cover 
the interval between the 16'^*' and the 84**' percentile. 


TABLE 1. Input parameters of the model 

a = 1.3 ±0.7 K"* 


A'a = 240 ± 90 mN m“* 

L = 9.8 nm 

AH° = — 54 ± 5 kcal moU* 
AS° = -154 ± 13 cal moU* K’* 


III. RESULTS AND DISCUSSION 
A. Qualitative observations 

In Fig. 3 we can visually compare confocal images of 
a DNA-functionalised GUV adhering to a SBL (a) and 
a non-adhering GUV on a passivated glass surface (b). 
Both GUVs and SBL are stained with fluorescent lipids. 
It is clear from both the 3D reconstruction and the 
vertical cross section that the adhering GUV takes the 
shape of a truncated sphere, with a flat and circular 
contact region. This evidence confirms the assumption 
that, at low enough temperature, DNA-mediated ad¬ 
hesion is strong enough to guarantee the dominance 
of stretching over bending and other contributions to 
the deformation energy. The fluorescence intensity 
measured within the adhesion patches is almost exactly 
equal to twice the value measured on the SBL outside 
the adhesion region (1.95 times for the vesicle shown 
in Fig. 3a). This evidence confirms the presence of two 
lipid bilayers in close contact within the adhesion area 
and excludes the possibility of DNA-mediated fusion of 
the two membranes. [52] 

The non-adhering GUV displayed in Fig. 3b does not 















:ontact angle, 6 (deg) 

O) 03 

O O O 

a 

1 

a:^ 10 
(/> 

1 ^ 
1 ^ 



20 40 60 

CL 4 

20 40 60 

^ 0.2 

d 

^ d 0.2 

e 

^ 0 
g-' -0.2 

^ -0.4 

O • 

• 

• 

•< 

20 40 60 

^ Q. 0 

cc 

g' -0.2 

S Q. 

E -0.4 

20 40 60 


Temperature, T (°C) Temperature, T (°C) 



Temperature, T (°C) 



FIG. 4. Experimental and theoretical temperature-dependent 
vesicle adhesion for samples with intermediate DNA coating 
density Pdna = 390 ± 90 fim~^ (see experimental section). 
In the top row we demonstrate the temperature-dependence 
of the contact angle 9 (a), adhesion-patch radius Rp (b), and 
vesicle radius R (c) for a typical adhering GUV. Points in¬ 
dicate experimental data, solid lines mark theoretical predic¬ 
tions, with errorbars visualised as grey-shaded regions. Fit¬ 
ting parameter Tq = —40°C. In the bottom row we sum¬ 
marise the results for 5 vesicles. Points represent the relative 
deviation of experimental data from theoretical predictions 

(^exp _^th)/^th, ^ ^ 

shaded regions mark the uncertainty interval of the theory. 
Model parameters are reported in Table I. 


FIG. 5. Experimental and theoretical temperature-dependent 
vesicle adhesion for samples with low DNA coating density 
Pdna = 39 ± 9 pm“^ (see experimental section). In the top 
row we demonstrate the temperature-dependence of the con¬ 
tact angle 6 (a), adhesion-patch radius Rp (b), and vesicle 
radius R (c) for a typical adhering GUV. Points indicate ex¬ 
perimental data, solid lines mark theoretical predictions, with 
errorbars visualised as grey-shaded regions. Fitting parame¬ 
ter To = —20°C. In the bottom row we summarise the results 
for 5 vesicles. Points represent the relative deviation of exper¬ 
imental data from theoretical predictions 
for X — 9 (d), Rp (e), and R (f). Grey-shaded regions mark 
the uncertainty interval of the theory. Model parameters are 
reported in Table I. 


exhibit a flat adhesion patch, as clear from the vertical 
cross-section. Note that for both the adhering and 
the non-adhering GUVs, the bottom part of the stacks 
appears brighter due to the z-dependent response of the 
instrument. 


B. Temperature-dependence of the geometrical 
observables 

In this section we discuss the temperature-dependence 
of the morphology of adhering GUVs. In Fig. 4a-c 
we show the experimentally determined contact angle 
0, adhesion-patch radius Rp, and vesicle radius R 
for a typical adhering GUV is a sample with DNA 
concentration equal to ~ ^90 ± 90 pm“^. The 

contact angle displays a non-monotonic behaviour as a 
function of temperature, with a positive slope at low T, 
followed by a sudden decay for T higher than ~ 30° C. 
The adhesion radius follows the same trend, as does 
the vesicle radius, which however displays much smaller 
relative variations. The solid lines in Fig. 4a-c represent 
theoretical predictions calculated using the input pa¬ 
rameters in Table. I, and using the neutral temperature 
To as a fitting parameter. Grey-shaded regions indicate 
propagated uncertainty in the theoretical predictions. 
The agreement between theory and experiments is 
quantitative at low temperatures. At high T, the theory 
fails to predict the drop in contact angle observed in 


experiments. This behaviour is expected since our 
theoretical description is valid in the limit of strong 
adhesion, where the attractive forces are sufficient to 
suppress bending contributions in the interaction energy. 
At high temperature the DNA, which in the present 
experiment features relatively short sticky ends, starts 
to melt, causing the loosening of the adhesive forces and 
a change in the GUV shape, detected as a shrinkage of 
the adhesion area. The temperature-dependence of the 
fraction of DNA bonds is quantified and discussed in the 
following sections. 

In Fig. 4d-f we show the relative deviations of the 
experimentally-determined morphological observ¬ 
ables from the theoretical predictions, defined as 
(Xe’^P _ for X = e, Ap, and A. The data, 

collected from 5 vesicles are consistent: the experimental 
data fall within theoretical errorbars at low T, deviating 
at higher temperature due to the failure of the strong 
adhesion assumption. 

In Fig. 5 we show experimental, and theoretically 
predicted morphological observables for the case of low 
DNA concentration, Pqna = 39 ± 9 /im“^. Similarly 
to the case of higher DNA concentration, experimental 
data are backed by theoretical predictions. In this case, 
however, the high-temperature deviation of the experi¬ 
ments from the theoretical predictions appears to be less 
evident, and shifted towards higher temperatures. This 
indicates the presence of an (however small) adhesive 
force hindering the partial detachment of the GUVs. 
We ascribe this behaviour to the effect of non-specific 
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membrane-membrane adhesion, e.g. dispersion attrac¬ 
tion, that for higher DNA coverage is suppressed by 
steric repulsion. 


C. DNA-melting 

We investigate the temperature-dependence of the 
fraction of formed DNA bonds via in-situ FRET mea¬ 
surements. Cyanine fluorophores, Cy3 (donor) and Cy5 
(acceptor), are connected to the 3' termini of a and a' 
sticky ends, as sketched in Fig. 1. FRET efficiency is 
described by 



where d is the distance between the fluorohpores and Rq 
is the Forster radius, equal to 5.4 nm for the case of Cy3- 
Cy5.[53] When sticky ends are bound to form a loop or 
a bridge, the distance between Cy3 and Cy5 is approxi¬ 
mately equal to the length of the hybridised sticky ends - 
2.7 nm, therefore we can assume a very high energy trans¬ 
fer efhciency for bound linkers. The average distance be¬ 
tween unbound linkers is sufficiently high to guarantee 
a comparatively very low transfer efficiency between un¬ 
paired tethers. Note that for FRET experiments the lipid 
membrane is not stained with Texas Red to avoid spuri¬ 
ous signal (energy transfer between Texas Red and Cy5). 
In Fig. 6a we show the confocal image of an adhesion 
patch (top), segmented to separate the actual adhesion 
area from the surrounding free SBL. This enables an ef¬ 
ficient characterisation of FRET efhciency in-situ. The 
emission spectra, collected within the patch and on the 
SBL while exciting the donor at 514 nm, are shown in 
Fig. 6c and d respectively. Colours from blue to red indi¬ 
cate low to high temperature. Spectra are normalised to 
the emission peak of Cy3, correctly found at ~ 568 nm. 
Note that for this experiment we used high DNA con¬ 
centration Pdna to strengthen the signal that otherwise 
would be too weak for a wavelength scan. As expected, 
the emission of the acceptor, peaked at ~ 665 nm, visi¬ 
bly drops at high temperature. We quantify this effect in 
Fig. 6d, where we plot the normalised acceptor emission 
intensity i = Icys/ (^Cys + ^Cy3),[54] where the /cy 5 /Cy 3 
are the peak-intensities estimated through a local Gaus¬ 
sian fit. Although qualitatively similar, the i curves mea¬ 
sured within and outside the adhesion patch exhibit some 
differences. For the case of free SBL, the emission ra¬ 
tio remains constant (~ 0.4) or slightly decreases upon 
heating, up to ~ 40° C, then it gradually drops down to 
~ 0.15. This decay is ascribed to the melting transition 
of DNA loops formed within the SBL. The FRET signal 
measured within the patch is higher at low temperatures 
(~ 0.6). This effect is probably due to the higher DNA 
density found within the patch at low temperature, which 


increases the probability of energy transfer between un¬ 
paired strands. Indeed we find that the overall fluores¬ 
cence intensity measured within the patch at T < 20° C 
is between 6 and 12 times higher than the intensity mea¬ 
sured on the SBL. The FRET signal measured within 
the patch is also found to increase upon heating, before 
suddenly decreasing at T ~ 45°C. The increase in FRET 
efhciency cannot be explained by an increase in DNA 
density within the patch, since the local DNA density de¬ 
creases as the adhesion area becomes larger upon heating. 
A possible explanation of this behaviour could be radia¬ 
tive cross-excitation between the two fluorophores, that 
becomes more efhcient as the adhesion patch gets less 
crowded upon heating. At high temperatures the FRET 
signal measured within the patch plateaus at ~ 0.15, in 
line with what we measure on the SBL. This confirms 
that at high enough temperature, when no bonds are 
formed, the DNA concentration is uniform across all the 
surfaces. 

The curves l{T) can be used to semi-quantitatively es¬ 
timate the temperature dependence of the overall frac¬ 
tion of DNA bonds. We Ht the low temperature plateaus 
(T < 35°C) in Fig. 6d with linear baselines B{T) and as¬ 
sume that l{T) plateaus to a constant value for T > 57°C. 
The fraction of formed DNA bonds is thereby estimated 
as 




l{T) {^t>57°C^ 
-®(^) ~ ('•T>57“C^ 


(17) 


A better estimate of (j){T) could be obtained by measur¬ 
ing t(r) up to higher temperatures, and fitting the high- 
temperature plateau with a second linear baseline. How¬ 
ever, temperatures higher than 65° C cannot be safely 
probed due to the risk of destabilisation of the dsDNA 
spacers. The experimental (('(T) data in Fig. 6e indi¬ 
cate that the DNA melting transition is relatively broad, 
spanning more than 30°C.[55] Moreover, the melting 
seems to occur at a higher temperature (by about 5°C) 
within the adhesion patch. 

The experimentally estimated fraction of DNA bonds can 
be compared with theoretical predictions. On the free 
SBL, only loops can form, and therefore ^(T) should 
be compared to the fraction of loops calculated 

accordingly to Eqs. 12, A17-A18. Within the adhesion 
patch we count contributions from bridges, loops formed 
on the GUV, and loops formed on the SBL. By assum¬ 
ing evenly distributed loops, the overall number of DNA 
bonds found within the patch is 


-Abound = N 



Xi +Xb 


-bcoxf^^Ap, 


(18) 


where x-t,/i are the fractions of loops and bridges on the 
GUV, given by Eqs. A27 and A28. The overall number 
of DNA tethers, bound and unbound, found within the 
patch is 


Ntot = N 



(1 - Xh) + Xh 


+ coAp, 


(19) 
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FIG. 6. In-situ FRET spectroscopy characterisation of the temperature-dependence of the fraction of DNA bonds, a, Confocal 
image of the adhesion patch of a GUV. At the bottom we show the same image segmented with our software to separate 
the adhesion patch (blue) from the free supported bilayer (green). Scale bar: 15 /rm. b, Fluorescence emission spectra 
measured within the adhesion patch for a typical vesicle (blue area in panel a). Golours from blue to red mark low to high 
temperatures in the interval 14.5 < T < 62.9°G. The curves are normalised to the emission peak of Gy3 (568nm). In panel 
c we show the emission spectra measured on the free SBL (green area in panel a), d, Relative intensity of the Gy5 emission 
peak (665nm) measured within (blue circles) and outside (green lozenges) the adhesion patch. The curves are relative to 3 
different vesicles. The amplitude of the Cy5 and Gy3 emission peaks is determined through a Gaussian fit of the 5 data points 
closest to the maximum, e, Experimental (symbols) and theoretically predicted (solid lines) fraction of formed DNA bonds 
within (blue circles) and outside (green lozenges) the adhesion patch. Experimental data are extracted from the curves in 
panel d as described in the text. Theoretical curves are calculated using the parameters in Table I, DNA-coating density 
Pdna ~ 1200 ± 300 To = —20, and Ro = 10 ^m. Note that the value of To does not significantly affect these quantities. 

Blue and Green shaded regions indicate propagated uncertainties of the blue and green solid lines. Gyan shaded region marks 
their overlap. 


where we assume that also unbound DNA is evenly 
distributed across the surfaces. The theoretically 
predicted fraction of bonds within the patch is thus 
.^boundZ-^tot- In Fig- 6e we compare theoretical ^(T) 
with experimental estimates. Since the choice of the 
neutral temperature Tg and radius i?o do not noticeably 
affect the melting curves, theoretical predictions are 
calculated using a fixed Tq = —20°C and i?o = 10 /rm, 
with no fitting parameters. Our model captures the 
width of the DNA transition as well as the difference 
in melting temperature between the patch and the free 
SBL. However the theory underestimates the average 
melting point by 5 — 10°C. This deviation is, at least 
partially, ascribable to the attractive effect of Cy3 and 
Cy5 molecules on the stability of DNA. For duplexes 
labeled with either of the molecules, the stabilisation 
has been quantified in a positive melting-temperature 
shift of 1.4 — 1.5°C.[56] The presence of both molecules 
is expected to cause a greater shift. Another explanation 
could be an underestimation of the DNA concentration 
Pqna- discussed below, this hypothesis is consistent 
with tension measurements. The impossibility of prob¬ 
ing the high-temperature baselines could also play a role. 


D. Membrane tension 


The temperature dependence of the membrane tension 
is measured by flickering analysis of the equatorial cross 
sections of GUVs. The tension is extracted by fitting the 
power spectra of the thermal fluctuations, determined as 
explained in the experimental section, with the function 


{\h\q.)\) 


ksT 

2Ca 


1 

Qx 



( 20 ) 


where C is the contour length of the equatorial cross- 
section of the GUV, a is the membrane tension, k the 
bending modulus, and qx the wave vector evaluated along 
the contour. Equation 20 is derived from the original 
work of Helfrich, [33] describing the fluctuations of an in¬ 
finite 2D membrane, and corrected to account for the 
fact that, by imaging an equatorial cross-sections, only 
modes propagating along the horizontal direction should 
be considered. [34] Of the discrete set of wave vectors 
qxin) = 27rn/£, modes with n < 6 are excluded from 
the analysis. Mode n = 0 and n = 1 correspond to size 
changes and translations of the GUV. Modes with n > 2 
describe thermal fluctuations. However, Eq. 20 is de¬ 
rived for a planar membrane, and should not be used to 
describe modes with n < 6, which are influenced by the 













11 




With the present technique we cannot access the tension 
of vesicles adhering to SBL for the case of higher DNA 
concentrations. Indeed, for values of a in the grey-shaded 
region on Fig. 7a, the relevant portions of the fluctuation 
power spectra are masked by experimental noise deriving 
from the finite resolution of the contour-tracking proce¬ 
dure. 

The membrane tension can be evaluated within the 
framework of our model. At equilibrium, the derivative 
of the interaction energy in Eq. 2, taken with respect to 
the GUV area, is 


dU _ dUuNA 

M dA 


= 0 , 


( 21 ) 


where we used 


FIG. 7. Experimental and theoretical temperature-dependent 
membrane tension, a. Membrane tension cr measured from 
flickering experiments. Blue circles indicate adhering GUVs 
with low DNA coverage (Pdna = 39±9 /rm“^). Red lozengges 
indicate non-adhering GUVs on a BSA-coated glass surface. 
Grey-shaded regions mark regimes in which membrane ten¬ 
sion is currently inaccessible to our technique. The blue- 
shaded band indicates the theoretical prediction for a cal¬ 
culated with the fitting parameter —60 < To < 0°G and low 
DNA coverage. Blue lines mark the corresponding errorbars. 
The green-shaded band indicate the theoretical prediction for 
a for higher DNA coverage (pdna = 390 ± 90 Green 

lines mark the corresponding errorbars. Model parameters 
are reported in Table I. [30] b. Fluctuation-amplitude spec¬ 
tra for adhering GUVs with low DNA coverage and various 
temperatures. Symbols indicate experimental data and solid 
lines indicate fits according to Eq. 20. From top to bottom 
the fitted values of the membrane tension are cr = 4.1 ± 0.5, 
2.42 ± 0.09, 1.88 ± 0.09 x 10“^ Nm. 


spherical geometry of the GUV. [34] For our analysis we 
fit the spectra for modes 6 < n < 16. At higher q we 
approach the resolution limits of the current method. 

In Fig. 7a blue circles mark the tension measured as a 
function of temperature for adhering GUVs. In Fig. 7b 
we show examples of power spectra fitted by Eq. 20. The 
tension typically lies in the interval 2 x 10“^ — 2 x 10“® 
N m“^, with clear variations between different GUVs. In 
the tested range, the tension consistently displays a weak 
dependence on temperature changes. 

For comparison, the membrane tension is measured on 
non-adhering GUVs supported by a passivated glass sub¬ 
strate. The values of a measured for non-adhering GUVs 
(red lozenges in Fig. 7a) are signihcantly lower than those 
measured on adhering vesicles, falling within the range 
10“® — 10“® N m“^, and being clustered around 10“® 
N m“^. Furthermore, membrane tension of non-adhering 
GUVs typically displays a strong decrease upon increas¬ 
ing temperature. The large variability observed in the 
tension of adhering and, in particular, of free GUVs, is 
ascribed to the polydispersity of electroformed samples, 
which produces vesicle populations with very different 
excess areas (Tq). 


_ ^Ustretching 

^ 5A 


= K, 


A-i 

A 


( 22 ) 


Equation 21 suggests that by measuring a we can di¬ 
rectly probe the DNA-mediated forces. The blue-shaded 
region in Fig. 7a marks the model predictions for u, 
calculated using the parameters Table I, = 39 ± 9 

/im“^, and values of the neutral temperature covering 
the experimentally observed range (—60 < Tq < 0°G). 
The size of the vesicles does not impact the predictions 
of (T, therefore we fix i?o = 10 pm. Solid blue lines mark 
the uncertainty interval propagating from the errorbars 
of the model parameters (Table I). With no fitting 
parameters, we observe a semi-quantitative agreement 
between theory and experiments. In particular, the 
theory predicts the weak temperature-dependence of a 
observed in the experiments. The slight underestimation 
of the theoretical predictions as compared to the mea¬ 
sured tension could be explained by an underestimation 
of the DNA-coating density. This evidence is consistent 
with the underestimation of the DNA melting tempera¬ 
ture discussed above. 

The green-dashed region in Fig. 7a indicates the theoret¬ 
ical prediction calculated using = ^90 ± 90 pm“^. 
The predicted tension falls within the non-accessible 
region. 


IV. CONCLUSION 

In this article we experimentally investigate 
temperature-dependent adhesion of Giant-Unilamellar- 
Vesicles on supported lipid bilayers mediated by mobile 
DNA linkers. The simple geometry of the problem allows 
for an accurate characterisation of the morphology of 
adhering GUVs and the temperature dependent fraction 
of bound DNA tethers by means of confocal microscopy. 
For the first time to our knowledge, we quantify the 
temperature-dependent membrane tension induced by 
DNA bonds by analysing the thermal fluctuations of the 
GUVs imaged across their equatorial plane. 
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The experimental results are compared to theoretical 
predictions from our recently developed model, [30] which 
we here extend to the case of vesicle-plane adhesion. The 
model takes into account both the elastic deformation 
of the GUV and the statistical-mechanical details of the 
DNA-mediated interactions. 

For sufficiently high DNA coverage, the adhesion contact 
angle exhibits a re-entrant temperature dependence. 
Upon heating from low temperature the contact angle 
increases, reaching a maximum at T ~ 30 — 40°C. 
Upon further temperature increase, the contact angle 
drops. The re-entrance is less pronounced or absent for 
lower DNA coverage. With a single fitting parameter, 
the model is capable of quantitatively predicting the 
low temperature regime and ascribes the increase in 
contact angle to the interplay between the temeperature- 
dependent excess area of the GUV and the entropic 
coupling between the hybridisation free-energy of the 
mobile tethers and the adhesion area. The theory is 
developed in the limit of strong adhesion, therefore it 
fails to predict the re-entrant behaviour of the adhesion 
area, caused by the weakening of the DNA bonds. 
The less-pronounced re-entrance observed for low DNA 
concentrations is ascribable to non-specific adhesive 
interactions that kick-in at high temperature, and are 
suppressed by steric repulsion for samples with high 
DNA coverage. 

The melting of DNA bonds is investigated in-situ by 
FRET measurements. We observe a broad melting tran¬ 
sition and find that bonds formed within the GUV-plane 
adhesion patch are more stable than in-plane bonds 
formed on free bilayers. With no fitting parameters our 
model can semi-quantitatively reproduce these features, 
although an underestimation of the melting temperature 
is observed. 

Membrane tension measurements performed on adhering 
GUVs demonstrate a weak temperature dependence. In 
a similar range of temperatures, non-adhering GUVs ex¬ 
hibit significantly lower tension, rapidly decreasing upon 
heating. The differences in magnitude and trend demon¬ 
strates the role played by DNA in mediating membrane 
adhesion. Experimental results are in semi-quantitative 
agreement with theoretical prediction, which further 
demonstrates the accuracy of the model used to describe 
hybridisation free-energy of tethered mobile linkers, and 
in particular the translational-entropic contributions 
that couples it to the adhesion area. 


Our experimental observations and the agreement with 
the theoretical predictions help to clarify the complex 
mechanisms controlling adhesion of soft units mediated 
by multiple linkers. Besides the fundamental interest for 
the still poorly understood physics of multivalent inter¬ 
actions, our findings can help the design of functional, re¬ 
sponsive, tissue-like materials with promised applications 
in biosensing, encapsulation-release mechanisms, and fil¬ 
tration. Finally, the conclusions drawn for our model 
system can be adapted to the quantitative description 


of cell adhesion and spreading on solid substrates, with 
possible biomedical implications in prosthetics and scaf¬ 
foldings for tissue regeneration. [57-61] 
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Appendix A: Details on model derivation 


1. Geometrical expressions 


In the limit of strong adhesion the GUV will take 
the shape of a truncated sphere with contact angle 9 
(Fig. 1), which we take as the independent variable of our 
model. [41] In this simple geometry, the contact (patch) 
area, total area, and volume of the GUV are respectively 

Ap = ttR^ sin^ 6 (Al) 

A = ttB? (1 -I- cos 9) (3 — cos 9) (A2) 

V = (1 -I- cos 9)^ (2 — cos 9). (A3) 

In the limit of water-impermeable membranes, the inter¬ 
nal volume of the GUVs can be taken as a constant 


V=\^Rl (A4) 

where we introduce a reference radius i?o- By using 
Eqs. A3 and A4 we obtain 


R — Rq 


4 

(1 -I- cosd)^ (2 


1/3 


COS0) 


(A5) 


which can be inserted into Eqs. A1 and A2 to make the 
0-dependence of Ap and A explicit. 

Let us now recall the definition of reduced volume of a 
vesicle [41, 62] 


3V 



By combining Eq. A6 with expression for the 
temperature-dependent unstretched area in Eq. 4, we ob¬ 
tain 


v=[l + a{T-To)]-^/\ (A7) 
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The reference temperature Tq is therefore defined as 
the temperature at which a GUV has reduced volume 
equal to 1. For T < Tq, when v > 1, an isolated vesicle 
resembles a turgid sphere, with non-zero membrane 
tension whereas for T > Tq, u < 1, it assumes the 
features of a “floppy” balloon, with excess area. At 
T = Tq an isolated vesicle is a perfect sphere with 
zero-membrane tension and radius equal to the reference 
radius i?o- 

By combining Eqs. A2, A5 and 4 we obtain an explicit, 
0-dependent expression for the stretching energy in 
Eq. 3. Note that Eq. A5 has been derived under a 
constant-volume assumption (Eq. A4). Alternatively, an 
equivalent relation can be derived for water permeable 
- solute impermeable - GUVs, in which the volume is 
set by the balance between the osmotic pressure drop 
across the membrane and the Laplace pressure. [41] In 
relevant experimental conditions the two assumptions 
lead to very similar results. [30] 


2. Free energy for bridge and loop formation. 

Eor the case of mobile linkers, the configurational en- 
tropic contribution to the bridge/loop formation free en¬ 
ergy (Eq. 5) can be split into a rotational and 

translational contribution 

= AA™* -b AS'‘7/". (A8) 

The rotational contribution takes the same expression for 
loop and bridge formation, and encodes for the reduction 
of configurational entropy following the hybridisation of 
two rigid tethers with fixed grafting sites [22, 30] 


Eqs. A8-A11 with Eq. 5, we obtain the hybridisation 
free energy of bridge formation on the loops formation 
on the GUV in Eq. 6. 

3. Fraction of loops and free tethers on the SBL. 

We now focus on the description of the tethers an¬ 
chored to the SBL, and calculate the equilibrium fraction 
of formed loops in the absence of an adhering GUV. This 
information is needed for the calculation of the GUV- 
SBL adhesive interaction as well as for a direct compar¬ 
ison with experimental data. 

Let us consider a finite portion of the SBL of area E, 
containing two populations of N linkers with a and a' 
sticky ends. Following Eq. 6, the free energy for loop (£) 
formation on the SBL can be written as 

= AG° - keT log . (A12) 

LpoAS 

By indicating as the number of loops within the SBL, 
and taking into account combinatorics, we can write the 
partition function of this systems as[30] 

^ = E (-mAGp^) , (A13) 

Ne ^ 

which can be rearranged as 

Z = (A14) 

Ne 

We now consider the limit of an infinite SBL with a con¬ 
stant DNA surface density cq, i.e. we take A, E —)• oo, 
with Co = A/E = pdna/ 2. By using the Stirling approx¬ 
imation we obtain 


where po = IM is a reference concentration. The trans¬ 
lational contribution encodes for the lateral confinement 
following the binding of two mobile tethers. Eor the case 
of loops, upon binding, two tethers initially capable of 
exploring the entire GUV surface area A, are confined to 
within a region ~ from each other [30] 

A5/-- = fcBlog =felog ^ . (AlO) 

For the case of bridge formation, we consider a free linker 
on the SBL, initially located within the contact area Ap, 
and a second linker on the GUV, which is free to ex¬ 
plore the entire surface area A. Upon binding, the area 
available to the pair is reduced to 47rL^Ap, resulting in 



S{ce) = const • E[-Cf logc^ - 2(co - Ci) log(co - Ci) 

— CiPAG^ — Ci log(poL) — Ci + const], (A15) 

where we define the density of loops as Ci = NijTi. 
Within the saddle-point approximation, [30] the sum in 
Eq. A14 is dominated by the stationary point of S 

g=0. (A16) 

By solving Eq. A16 we obtain the expression in Eq. 12 
in the text, where we introduce the concentration of free 
tethers on the SBL Cf = cq — c^. Note that with Eq. 12 
we recover a simple mass-balance relation between loops 
and free tethers on the SBL, which ultimately results in 



(A17) 

(A18) 


We notice that, contrary to the case of two adher- The fraction of loops formed within the bilayer is = 

ing GUVs, [30] here AS'J]''®'”® = AS')"'“®. By combining c^/cq. 
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4. Combinatorial effects. 


pressed as 


Given the expressions for the hybridisation free-energy 
of a single bridge/loop (Eq. 6), a combinatorial approach 
is required to compute the overall DNA-mediated inter¬ 
action energy. 

Following the derivation carried out to describe loop for¬ 
mation on the SBL, we indicate the total number of teth¬ 
ers with a {a') sticky ends on the GUV as N, and define 
N( as the number of those tethers linked in loops. We 
indicate as iVbi, with i = 1,2, the number of tethers 
forming bridges with those on the SBL, with the index 
i referring to a and a' sticky ends. We label as Na the 
number of free tethers on the SBL located within the ad¬ 
hesion patch, where the index i = 1,2 now refers to a' 
and a sticky ends. The partition function of the system 
of linkers can be written as 

z{Nfi, Nfi, N) = ^ [^Nfi,N{2,NiNhi,Nbi,N() 

Ni,„Ni 

exp [-/3{Nhi + Nh 2 + Ne)AG ]], 


A = l3AG*Xi + xi (log a;^ -I-1) -I- 

^ |^/3AG*a;b* + xu (logiV - logVf - l) + 

2^1,2 

Xbi (logXbi + 1) + {XH - Xhi) log {xii - Xhi) 


+ {1- xe- Xbi) log {1- xe- Xhi) 


(A24) 


Note that in Eq. A24 we re-defined the hybridisation free 
energy for bridge and loop formation as 

AG* = AG- ksT log N. (A25) 


For typical experimental conditions, the attractive com¬ 
binatorial term —k^T log N in Eq. A25 can be estimated 
in « — lOfceT.pO] Within the saddle-point approxima¬ 
tion, the sum in Eq. A23 is dominated by the stationary 
point of A 

DA 

-— = 0 with ?/=bl, b2, 1, fl, f2. (A26) 

OXy 


where the number of possible configurations for a given 
Ni and iVbi is 


^Nn,N{2,N{A!hl,A[\2i,Ni) — 



iVb^! 


(A19) 


From the saddle-point equations Eq. A26 we obtain 
Eqs. 8 and 9 in the text, where we find a;bi = Xh 2 = a^b- 
By solving Eqs. 8 and 9 we obtain 


9b 


Xb = 


+ 2<7b + 4(jf^ -b 1 — (7b ~ 1 
2(7t 


(A27) 


To account for strand-concentration fluctuations within 
the adhesion patch, we need to consider that Na is 
Poisson-distributed around its average value iVf 

P{Nu,Ni) = exp[-Af]^. (A20) 

iVf^! 

Using Eq. A17, and recalling that Cf = cq — Q is the 
concentration of free tethers within the SBL, we find 

iVf(0) =cfAp(0), (A21) 


9b + 29b + 2(7£ -b 1 — (gb + 1) ^9^ + 2gb + 4g^ -b 1 

Xt = - 7, -• 

29t 

(A28) 

Note that for simplicity the fraction of bridges and loops 
are indicated as Xb/^. The saddle point equations for Xf 
(xfp = Xf_ 2 ) read Xf — Xb = N[/N, which confirms that 
the density of the free tethers in in the patch region is 
equal to that of the reservoir, as expected. 

By inserting the saddle-point solutions for x^, Xb, and Xf 
in Eqs. A23 and A24 we can calculate the free energy 
Uhyb (Eq. 7).[19] 


where we highlighted the strong dependence on the con¬ 
tact angle 9. Using Eqs. A19 and A20 we write the full 
partition function as 


ZNn,NniNt,N) =Y,PiNii,Nt)P{N!2,Nf)z{Na,Nn,N). 


Nt 


Equation A22 can be rearranged as 


(A22) 


5. Reference energy. 

The reference energy Uq in Eq. 2 is calculated for iso¬ 
lated GUV and SBL and can be written as 

Uo = C/g^tretching ^ (^29) 

The stretching term is [30] 


E exp[—AA]. (A23) 

-^fl ,^hl ,^h2 


xrstretching 

^0 


0 if T > To 


(A30) 


By defining the fractions Xy = Ny/N {y =bl, b2, 1, fl. Note that the stretching contribution is only present for 
f2), and using the Stirling approximation, A can be ex- pre-stretched vesicles, i.e. if the reduced volume is v > 1 
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(i.e. T < To). [41] The DNA contribution is calculated for 
a GUV of area equal to the unstretched area A, in which 
only loops can form. By following the steps outlined in 
section A 4 and in ref. [30] we calculate the fraction of 
tethers involved in loops 


0 + 1 — V 4g^ + 1 

= - ’ 


(A31) 


where qi is given by Eq. 10. The DNA part of the refer¬ 
ence energy is 


^ jv/ceT x'l + 21og (l - - 2^ - 2log j 

(A32) 

where N{ = C{Ap is the number of free tethers present 
within area Ap on the SBL. Ap is the zero-stretching ad¬ 
hesion area, which the GUV-SBL system would form for 
negligibly small attractive forces when T > Tq, as de¬ 
rived in ref. [30] 

Note that Uq does not depend on the contact angle 9 
therefore its form does not influence the equilibrium fea¬ 
tures of the system. 
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